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Abstract. We use three samples (3CRR, 6CE and 6C*) to investigate 
the radio luminosity function (RLF) for the 'most powerful' low-frequency 
selected radio sources. We find that the data are well fitted by a model 
with a constant co-moving space density at high redshift as well as by one 
with a declining co-moving space density above some particular redshift. 
This behaviour is very similar to that inferred for steep-spectrum radio 
quasars by Willott et al (1998) in line with the expectations of Unified 
Schemes. We conclude that there is as yet no evidence for a 'redshift cut- 
off' in the co-moving space densities of powerful classical double radio 
sources, and rule out a cut-off at z ;$ 2.5. 



1. Introduction 

It is well-known that the (co-moving) space densities of the rarest, most powerful 
quasars and radio galaxies were much higher at epochs corresponding to z ~ 2 
than they are now (Longair 1966). The behaviour of the space density beyond 
these redshifts is the subject of this paper. Dunlop & Peacock (1990) found 
evidence for a 'redshift cut-off' (a decline in the co-moving space density) in 
the distribution of flat-spectrum radio sources over the redshift range 2 — 4. 
Through failing to find any flat-spectrum radio quasars at z > 5 in a large (~ 40 
per cent of the sky) survey, Shaver et al (1996, hereafter SH96) argued for an 
order-of-magnitude drop in space density between z ~ 2.5 and z ~ 6, for this 
class of object. As emphasised by SH96, the crucial advantage of any radio- 
selected survey is that with sufficient optical follow-up, it can be made free of 
optical selection effects, such as increasing dust obscuration at high redshift. It 
is chiefly for this reason that the work of SH96 provides the most convincing 
evidence to date for the existence of an intrinsic decline in the co-moving space 
density of any galaxy class at very high redshift. 
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Figure 1. L151 — z plane for the three samples used in our analysis. 
3CRR (diamonds; Laing, Riley & Longair 1983), 6CE (crosses; Eales 
et al 1997), 6C* (stars; Blundell et al 1998; Jarvis et al in prep.) and 
Cygnus-A (filled triangle). The filled circles represent the three 6C* 
sources with upper redshift limits and the associated arrows show the 
paths of the possible redshifts and luminosities for these sources. The 
area between the horizontal lines is the region which contains the 'most 
powerful' sources according to our definition. The dark grey region 
shows the approximate L151 — z plane coverage of the 3CRR sample 
(10.9 Jy < S178 < 80.0 Jy, spectral index, a = 0.5) and the light grey 
region the 6C* sample (0.96 Jy < S151 < 2.0 Jy, with a = 0.98 and 
a = 1.5 for the lower and upper fluxes, respectively). Note the area 
between the two shaded regions contains no sources, this is the area 
which corresponds to the absence of a flux-limited sample between the 
6CE (5151 < 3.93 Jy) and 3CRR (S 178 > 10.9 Jy) samples. 



2. Modelling the RLF 

We adopt a parameterisation of the RLF which is separable in 151-MHz lumi- 
nosity L151 and redshift z with a single power-law in L151 of the form [Lj X )~ /3 . 
We consider two cosmologies f^M = 1> = (cosmology I) and Qm = 0.3, 
f^A = 0.7 (cosmology II). Model A parameterises the redshift distribution as a 
single power-law of the form (1 + z) n . For model B the redshift distribution 
is parameterised as a Gaussian, giving an overall expression for the co-moving 
space density of p = p {L/ L )~" exp —{(z — z Q )/y/2 z\} 2 where p and [3 are the 
normalising term and power-law exponent respectively, z Q is the Gaussian peak 
redshift and z\ is the characteristic width of the Gaussian. Model C is described 
by the same model up to z Q beyond which it becomes constant. 
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3. Results and Discussion 



For sources in the top-decade in luminosity of the L151 — z plane (Fig. 1) our 
parametric fitting and likelihood analysis of model radio luminosity functions 
(Table 1) show that the data are inconsistent with a (1 + z) n power-law in 
redshift (Model A), but are well fitted by both models B and C. These models 
are shown in Fig. 2 in the form of a log N / log S plot. We conclude that although 
the relative likelihood for model B is 2.5 times larger than for model C, this is not 
statistically significant enough to distinguish between the two models with any 
confidence. This uncertainty is further compounded by the effects of assuming a 
mean spectral index in the model fitting. This result is in very close agreement 
with the RLFs for radio loud quasars modelled by Willott et al (1998) and 
various studies of AGN at optical (Irwin et al 1991) and X-ray (Hasinger et al 
1998) wavelengths. 

This is in apparent contradiction to the findings of SH96 for flat-spectrum 
quasars. If the relationship between the flat- and steep-spectrum populations is 
as described by unification models of AGN then we might expect to see similar 
evolution in the two populations. Thus to determine the co-moving space density 
of radio sources at high-redshift, an understanding of the spectral index trends, 
K— corrections and associated selection effects must first be achieved. 

Fig. 2 also illustrates the contribution of powerful sources at high redshift 
to the total source count in a low-frequency survey. We see that even for the 
no cut-off model (Model C) the fractional contribution is very small. This may 
render the location of the redshift cut-off virtually impossible to determine until 
the selection effects associated with radio surveys are fully understood. 
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Table 1. Best-fit parameters for the model RLFs, described in the 
text. Pks is the 2-D Kolmogorov-Smirnov probability (a value above 
0.2 signifies a reasonable fit to the data) and C-ji is the likelihood rela- 
tive to model B. 

Errors for model B (cosmology I): Alog 10 (p o ) = 0.17, A/3 = 0.2, 
Az = 0.10, Azi = 0.17, for 68% confidence regions. 
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Figure 2. For cosmology II we show the areal density of the most 
powerful radio sources as a function of 151-MHz flux density S151 / Jy 
in comparison with the total source count. The 6CE source count 
(squares) is reproduced from Hales, Baldwin & Warner (1988) and 
the 7C source count (diamonds) from McGilchrist et al (1990). The 
3C source (stars) count was inferred from the revised 3CR sample of 
Laing, Riley & Longair (1983). The large open star represents the 6C* 
data (note that this is a lower-limit due to spectral index and angular 
size selection in this filtered sample); the open circle the 6CE sample; 
the inverted triangles the 3 bins representing the 3CRR sample. The 
solid horizontal lines show the S151 range of each sample. The lower 
panel shows the percentage of the total source count contributed by 
each model. Models A, B and C are represented by the dotted, dashed 
and solid lines respectively. 
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